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Introduction 

I N the hypersonic flight regime, the air-breathing supersonic 
combustion ramjet (scramjet) has been shown to be a viable 
propulsion system. The current designs of scramjet engines 
provide performance benefits only up to a Mach number of 
14.' Performance losses increase rapidly as the Mach number 
increases. To extend the applicability of scramjets beyond 
Mach 14, research is being conducted in the area of inlet and 
wave drag reduction, skin-friction and heat-transfer reduction, 
nozzle loss minimization, low-loss mixing, and combustion en- 
hancement. 

For high Mach number applications, hydrogen is the obvious 
fuel choice because of its high energy content per unit mass 
in comparison with conventional fuels. These flight conditions 
require engines to operate at supersonic internal velocities, 
high combustor temperatures, and low static pressures. The 
high static temperature condition enhances the production of 
radicals such as H and OH, and the low-pressure condition 
slows the reaction rates, particularly the recombination reac- 
tions. High-temperature and low-pressure constraints, in com- 
bination with a small residence time, result in a radical-rich 
exhaust gas mixture exiting the combustor. At high Mach num- 
ber conditions (due to low residence time), H and OH do not 
have enough time to recombine; thus, a significant amount of 
energy is lost as these high-energy free radicals are exhausted. 
Therefore, any enhancement of recombination reactions will 
aid in recovering a part of the energy that would otherwise be 
lost. 

One approach to enhancing the efficiency of hypersonic air- 
breathing propulsion .systems is the manipulation of the chem- 
istry by using additives, sensitizers, or excited species in the 
combustor. Twarowski* showed that the presence of small 
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amounts of phosphine combustion products in a combustion 
mixture containing H and OH radicals increases the recombi- 
nation rates significantly at 1970 K. In that study, Twarowski" 
also assembled a mechanism for the reactions of phosphorus 
oxide and acid molecules, and estimated their reaction rates. 
Twarowski,^ in his follow-up study, showed that for fuel-rich 
combustion, the addition of 2— 3 wt% of phosphorus produced 
a specific impulse gain of 3 — 10 s for resp>ective initial reactant 
temperatures of 15(X)-1(XX) K. 

The objective of the present study is to conduct a flowfield 
analysis for a typical nozzle geometry fora NASP-type vehicle 
to assess the prospects for thrust enhancement in hypervelocity 
nozzles by substituting small amounts of phosphine for hydro- 
gen. A reaction sensitivity study of Twarowski’s^ recently de- 
fined phosphorus chemistry is also conducted to identify the 
critical rate-limiting reactions. 

Governing Equations and Flow Conditions 

A quasi-one-dimensional Euler code written specifically for 
reacting-flow problems is used to analyze the flowfield. The 
governing equations are solved numerically with a three-stage 
Runge-Kutta method and a third-order upwind-biased scheme 
based on Roe’s flux-splitting algorithm. The finite rate chem- 
ical reaction of phosphine, hydrogen, and air is modeled by a 
14-species 33-reaction model, which was obtained from Ref. 
2 and given in Table 1. The following flow conditions, which 
correspond to Mach 18 flight conditions, were chosen as the 
baseline conditions for the flow exiting the combustor at the 
nozzle throat: velocity = 5364 m/s, temperature = 3072 K, 
pressure = 101,325 N/m\ and fuel/air ratio (<^) = 1.85. 

The total mass flow rate of fuel is kept constant. When phos- 
phine is added to the fuel, it is substituted for hydrogen by 
mass. The total volume of fuel decreases as the percentage of 
phosphine increases because of its high density. A rectangular 
nozzle, representative for a 60-m single-stage-to-orbit vehicle, 
is considered in this study. The nozzle has an inflow area of 
0.1524 m^ and expansion ratio of 40. The area expands rapidly 
with an expansion angle of 24 deg for first 20% of nozzle 
length, and the expansion angle then reduces to 8.77 deg. 

Characteristic Time 

All chemical reactions have an associated characteri.stic time 
for which many definitions exist. A common measure is the 
time required for a major species involved in the reaction to 
reach a percentage of its equilibrium concentration. This def- 
inition requires the integration of a detailed chemical system 
from initial conditions to equilibrium. Because the primary 
motivation of this work is to determine chemical additions that 
decrease the characteristic reaction time of the hydrogen and 
air system, a precise measure of the characteristic time is 
needed to quantify the reaction system. 

The chemical Jacobian is defined by dwjdfj, where w, is the 
species production rate, and / is the mass fraction of species 
/. It correlates the effects of species concentrations and reaction 
rates, and includes the entire spectrum of rate information that 
ranges from the rate-limiting reactions to those that react at 
the fastest rate. Transformation of the chemical Jacobian into 
diagonal form with the appropriate similarity transform yields 
the eigenvalues of the chemical reaction. Each eigenvalue cor- 
responds to the characteristic time exhibited by a particular 
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Table 1 Reaction set and rate constants* 


No. 

Reaction*’ 

log, A 

n 

E 

h 

1 

H + H + A/^H. + Af 

17.81 

-1.0 

0 

— 

2 

O + O + A/^Oj + M 

17.00 

-1.0 

0 

— 

3 

H + O-fAf^OH + A/ 

16.42 

-0.6 

0 

— 

4 

H + OH + A^ ^ H 2 O + M 

21.92 

-2.0 

0 

— 

5 

H -f 0 , + A/ ^ HO 2 + M 

17.85 

-0.8 

0 

— 

6 

H -F 0 , ^ OH + 0 

17.20 

-0.927 

70.6 

— 

7 

OH 4 H. - H.O + H 

13.80 

0 

33.6 

— 

8 

OH 4 OH ^ d + H 2 O 

8.320 

1.4 

-1.7 

— 

9 

0 4 H 2 - OH 4 H 

14.27 

0 

57.0 

— 

10 

HO, 4 H -- OH 4 OH 

14.34 

0 

5.9 

— 

n 

HO 2 4 H -- H 2 4 0 . 

13.4 

0 

2.9 


12 

HO. 4 0 - O 2 4 OH 

13.3 

0 

0 

- — - 

13 

HO 2 4 OH ^ H.O 4 O 2 

13.3 

0 

0 



14 

0 . 4 PO 4 Af ^ PO, 4 A/ 

13.5 

— 


0 

15 

H 4 PO, 4 M - HOPO 4 M 

18.6 

— 

_ — ^ 

8.86C-04 

16 

H 4 PO 4 A/ HPO 4 M 

16.7 


— 

5.75C-04 

17 

0 4 PO 4 A/ - PO 2 4 M 

17.7 

— 


7.15C-04 

18 

0 4 PO. 4 A/ - PO, 4 Af 

18.2 

— 

— — 

9.23e-04 

19 

0 4 HOPO 4 A/ - HOPO. 4 Af 

18.2 

— 

— 

7.85e-04 

20 

OH 4 PO 4 M HOPO 4 M 

16.9 

— 

— 

4.80C-04 

21 

OH 4 PO 2 4 A/ - HOPO: 4 M 

20.5 

— 


1.19e-03 

22 

OH 4 HPO -- H 2 O 4 PO 

U.5 

0 

7.5 


23 

0 4 HOPO - OH 4 PO. 

13.5 

0 

8.3 

— 

24 

H 4 HOPO H.O 4 PO 

13.5 

0 

24.4 

— 

25 

H 4 HOPO: ^ H 2 O 4 PO. 

13.5 

0 

0.2 

— 

26 

H.O 4 PO, -- OH 4 HOPO. 

11.5 

0 

0 

— 

27 

OH 4 PO ^ H 4 PO. 

11.5 

0 

23.3 


28 

0 , 4 PO ^ 0 4 PO: 

11.5 

0 

22.2 

— 

29 

0 4 HOPO: - 0 . 4 HOPO 

13.5 

0 

37.5 

— 

30 

H 4 HOPO ^ H. 4 PO. 

13.5 

0 

0.2 

— 

31 

OH 4 HOPO H.O 4 PO: 

1 1.5 

0 

0 


32 

OH 4 HOPO ^ H 4 HOPO: 

11.5 

0 

16.7 

— - 

33 

PO 4 HOPO: -- PO. 4 HOPO 

11.5 

0 

0 

— 


“Units are kilojoules per mole, seconds, moles, cubic centimeters, and degrees Kelvin. 

^Third-body efficiencies for all thermolecular reactions are 2.5 for M = H 2 , 16.0 for H;0, and l.O for 
all other M. Kf = AT" exp(-t7/?T) for reactions I - 1.1 and 22-33, and 410 for reactions 14 21. 


reaction within the system. The largest and smallest eigenval- 
ues correspond to the shortest and longest characteristic times, 
respectively. The overall reaction system is controlled by the 
rate-limiting reaction, which has the longest characteristic 
time. 

Thus, to a first-order approximation, the smallest eigenvalue 
(An, in) of the chemical Jacobian provides a measure of the over- 
all limiting reaction rate of the chemical system, and the re- 
ciprocal of this eigenvalue determines the characteristic time 
tr of the system (r, = l/A^nn). In general, by perturbing any of 
the independent thermodynamic or species variables, and re- 
calculating the smallest eigenvalue of the chemical Jacobian, 
we can establish the sensitivity of the characteristic time to the 
perturbation. The relative importance of the reactants and the 
products can thus be identified, as well as the necessary degree 
of precision for the reaction constants. Finally, the influence 
of chemical additives can be quantified by correlating the 
change of characteristic time with incremental changes of ad- 
ditive concentration. 

Results and Discussion 

A sensitivity study of the reaction mechanism used in the 
study is performed to identify the critical reactions and an 
assessment of their reaction rate uncertainty on the accuracy 
of the prediction is made. The results indicate the following; 

Reaction 4: 

H + OH + M ^ H 2 O + M 

is the most sensitive reaction. A 10% perturbation changes the 
limiting (longest) characteristic time by 3.8%, which is the 
largest relative change among the 33 reactions. The second 
and third most sensitive reactions are, respectively, shown. 


Reaction 30: 

H + HOPO -- H. F PO. 

Reaction 15: 

H + PO. + A/ -- HOPO + M 

For these reactions, the characteristic time is changed by 2.5 
and 1.0%, respectively. The rate constants for these reaclit)ns 
must be well established for accurate nozzle performance pre- 
dictions. 

Figures 1 and 2 show the ratio of limiting characteristic 
times as a function of temperature and percentage of phosphine 
in the fuel for = 1.85 (fuel rich) and 4> ~ 0.6 (fuel lean) 
conditions. The ratio is the longest characteristic time 

without phosphine added to that with phosphine added. The 
efficiency of catalysis increases as the value of r, (,//, increases 
because of the decrease in t, for the rate limiting reaction. The 
ratio decreases as the temperature increases, and the abil- 
ity of phosphine to enhance the recombination decreases as the 
temperature approaches about 3000 K. This trend is in oppo- 
sition to the desired result in the scramjet application. An op- 
timum incremental amount of phosphine addition exists that 
gives the maximum decrease in r, . Figure 1 show's that for the 
typical high Mach number fuel-rich (<^>= 1.85) combustor-exit 
conditions, such as those considered in the present study, the 
ratio r, .,)//, is higher than fuel-lean conditions shown in Fig. 2. 
It shows that addition of phosphine w'ill be more effective at 
fuel-rich conditions as opposed to fuel-lean conditions. 

Flow computations were performed for several flow condi- 
tions. The baseline case was computed first, and then a para- 
metric study was conducted by systematically varying the in- 
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<^- 1.85 

r 

1 i4.ei«x 
D 13.0317 
C 13.03t3 
» ii.oato 
A 10.0244 

0 ». 03104 

• 0.0109 

7 7.01704 

4 4.01442 

9 0.01210 

4 4.00079 

> 3.00731 

2 2.00407 

1 1.00244 


Fig. 1 Ratio of longest characteristic time with and without 
phosphine as function of temperature, and the percentage of phos- 
phine in fuel reaction rates for tff = 1.85. 



0 - 0.6 

r 3.24S0 

■ 2.14039 
D 2.07400 
C 1.00930 

■ 1.00300 

A 1.01030 
9 1.73200 

• 1.44737 

7 1.94107 

C 1.47437 
S 1.39007 
4 1.30934 

3 1.21006 

2 1.13434 

1 1.04005 


Fig. 2 Ratio of longest characteristic time with and without 
phosphine as function of temperature, and the percentage of phos- 
phine in fuel reaction rates for = 0.6. 



Fig. 3 Axial thrust with percentage of phosphine in fuel as pa- 
rameter for * 1.85. 



Fig. 4 Axial thrust with 5% phosphine in fuel with as a pa- 
rameter. 


flow conditions. Axial distribution of normalized thrust 
because of the addition of phosphine at 0 = 1.85 is shown in 
Fig. 3. The subscripts p and 0 refer to the thrust Th with and 
without phosphine addition, respectively. In the initial 20% of 
the nozzle section, the thrust increases very rapidly because of 
high flow expansion. As the amount of phosphine is increased, 
thrust production by the nozzle increases. It shows that 10% 
of phosphine addition enhances thrust by about 4.5%. This is 
a significant enhancement, given the sensitivity of thrust, to 
vehicle performance. 

Figure 4 shows the change in thrust caused by the 5% of 
phosphine addition, with (j> ranging from 0.5 to 3.0; moreover, 
it shows that thrust enhancement is a strong function of (f) and 
nozzle geometry. For a 10-m nozzle maximum, thrust is pro- 
duced for <^ = 2, but on the other hand, if the nozzle is greater 
than 20 m, = 1 produces maximum relative thrust. 

Conclusions 

A numerical study that used the quasi-one-dimensional Euler 
equations along with species conservation equations was con- 
ducted to study the possibility of thrust enhancement in hy- 
personic nozzles by substituting a small amount of phosphine 
for hydrogen to increase the recombination of radicals. 
The study showed that thrust production can be enhanced by 
phosphine addition. For nozzle inflow conditions correspond- 
ing to Mach 18 flight, the thrust can be enhanced by about 
4.5%. 
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